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Abstract. Spin-dependent optical potentials allow us to use microwave radiation 
to manipulate the motional state of trapped neutral atoms (Forster et al. 2009 
Phys. Rev. Lett. 103, 233001). Here, we discuss this method in greater detail, 
comparing it to the widely-employed Raman sideband coupling method. We provide a 
simplified model for sideband cooling in a spin-dependent potential, and we discuss it 
in terms of the generalized Lamb-Dicke parameter. Using a master equation formalism, 
we present a quantitative analysis of the cooling performance for our experiment, which 
can be generalized to other experimental settings. We additionally use microwave 
sideband transitions to engineer motional Fock states and coherent states, and we 
devise a technique for measuring the population distribution of the prepared states. 
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1. Introduction 

Motional state control of atomic particles is achieved by the absorption and emission 
cycles of a resonant or near resonant radiation, i.e., by light scattering typically at optical 
frequencies. For instance, laser Doppler cooling reduces the momentum of atoms or ions 
through multiple recoil processes [1]. Coherent momentum transfer can be performed 
with two-photon Raman processes [2] for applications in, e.g., atom interferometry [3]. 

The quantum state of the atomic particles is composed of the internal states, e.g., 
two spin states , \\)} for a two-level atom, and the external motional state. For 
free particles the simplest motional state is the momentum state \p). Trapped particles 
are instead characterized by vibrational eigenstates \n), which in the simplest case of a 
harmonic oscillator of frequency u; v ib have their energies equally spaced as hu) v ib(n+l/2). 

In free space, the momentum state of a particle, and consequently its kinetic 
energy, is changed by the momentum transfer \p) — > \p') in the absorption/emission 
cycle of an optical photon. While the momentum transfer picture also applies approx- 
imately for trapped particles when the energy separation between motional states is 
not spectroscopically resolved, recoil-free transitions become possible in the resolved- 
sideband regime (Mossbauer effect). While carrier transitions do not change the 
vibrational quantum state |n), the motional state can be controlled via sideband 
transitions \n) — > \n') (n' ^ n), for instance, in incoherent cooling processes \n) — > 
\n — 1) or in coherent manipulation of vibrational states [4]. With trapped ions or 
neutral atoms trapped in optical lattices, the resolved-sideband regime is typically 
realized by two-photon Raman transitions connecting two different hyperfme ground 
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states [5, 6, 7]. Alternatively, in spin-dependent potentials it becomes also possible to 
use microwave transitions, which also offer sufficient spectral resolution [8, 9, 10, 11, 12]. 

In the semi-classical picture shown in figure (1), an atomic transition exchanges 
either kinetic or potential energy with the motional degree of freedom of the atom. With 
the absorption of an optical photon, the kinetic energy is changed by the momentum 
kick from the photon, and quantum mechanically the process can be interpreted in terms 
of a displacement of the wavefunction in momentum space. With the absorption of a 
microwave photon, which carries a negligible momentum, the potential energy of the 
atom can be changed if the trapping potential of the two states are different; this allows 
an interpretation in terms of a wavefunction displaced in position space. 

2. Microwave induced motional sideband transitions 

2.1. Motional states in a state dependent lattice 

We consider a single atom with two spin states {|t)> trapped in a one-dimensional 
optical lattice. We will initially ignore the internal degree of freedom of the atom and 
take the Hamiltonian governing its motion in the trap as given by 



with Uq being the trap depth, k-^ = 27t/Al being the wavenumber of the two 
counter propagating laser fields creating the lattice, and x, p, the atom's position and 
momentum, respectively. 



Figure 1: In a semi-classical picture, an atomic transition can affect the motional state 
of an atom either (a) by a kinetic energy change caused by the momentum transfer 
from an optical photon of wavevector k opt (velocity selective transition [1]), or (b) by 
a potential energy change when the potentials of the two internal states are displaced 
in space by Ax (position selective transition). In the two cases, the motional energy is 
decreased when the detuning is set to the Doppler shift /c pt^at or the potential energy 
m u 2 Ax 2 / (2h) , respectively. 




(1) 



(a) 
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The motional eigenstates of an atom in such a potential are the well-known Bloch 
wavefunctions j^^), where n is the band index (n = for the first band) and k is the 
wavevector in the first Brillouin zone (BZ). In the limit of deep lattice potentials that we 
are considering here, the atoms remain localized for the time scales of the experiment 
and their spatial state is best described by the maximally localized Wannier state [13] 

M = ^=J2 e ~ lkrd \<*)- ( 2 ) 



N 

keBZ 

Here N is the lattice size, r the site index, and d = Al/2 the lattice spacing. In this 
deep lattice regime, we can safely view the vanishingly narrow energy bands e n (k) as 
the vibrational level energies e n of the corresponding Wannier state \n, r) at lattice site 
r; in the harmonic approximation we would have e n = hw^^n + 1/2). 

The Wannier states form an orthonormal basis set such that the overlaps between 
two different states yield (n, r | n',r') = S nin 'S r y. This means that the interaction of the 
atomic spin with a microwave field will fail to induce motional sideband transitions, 
\n, r) -H- \n',r'), because of the nearly negligible momentum carried by microwave 
photons, five orders of magnitude smaller than that by optical photons. This restriction 
can be lifted if the atom experiences a different trapping potential depending on 
its internal spin state as the corresponding motional eigenstates are then no longer 
orthogonal [14, 15]. A simple relative spatial shift of the potentials trapping each internal 
state induces such a difference. A shift by a distance Ax is accounted for by the position 
space shift operator T& x = exp(—ipAx/h), see figure (2a). The overlap between the 
two Wannier states then becomes 

Tax 



n,r 



n,r> Ax). (3) 



The resulting overlap integral, — 1 < I™£ (Ax) < 1, is hence a known function of Ax, 
see figure (2b). It is analogous to the Franck-Condon factor from molecular physics and 
it determines the strength of the transitions coupling different vibrational levels [16]. 

One way to realize the shift operator T& x is by two overlapped lattices which trap 
each spin state separately and can be independently shifted in the longitudinal direction 
as shown in figure (2). The trapping potential thus becomes dependent on the spin state 
s = { j, 1} and the shift distance Ax = x° — a;?, 

# ext = f^+ E ^cos 2 [fc L (x-a;°)]®| S )( S | (4) 

with x° s being the position of the lattice trapping the state \s). The total transition 
matrix element for two spin states coupled by an interaction Hamiltonian Hj, with a 
free-atom bare Rabi frequency Q Q , is then given by 

hn<;'(Ax)/2 = (s\ n\ r' \f Ax ® ifJ s, n, r) = 1$' (Ax) x m /2. (5) 

The Franck-Condon factors J^jT (Ax) can be explicitly evaluated using equations (2) 
and (3). We first rewrite equation (2) using Bloch's theorem, 

W n , r (x) = ^=J2 E e~* krd e^ a n>q (k) \k) . (6) 
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Figure 2: (a) The coupling strength of a sideband transition in a spin-dependent lattice 
is the bare spin state coupling VIq multiplied by the overlap between the two involved 
vibrational states, the Franck-Condon factor, which is controlled by the relative shift Ax 
between the two lattices. rj x is the spatial Lamb-Dicke parameter defined in section 2.1 
and later in 3.1. (b) Lattice shift dependence of the Franck-Condon factors for different 
transitions, denoted asn-m, calculated for typical experimental parameters (see text). 



with a n>q (k) being the Fourier coefficients of the Bloch functions and \k) the planewave 
state. These functions can be constructed using the periodic solutions of the Mathieu 
differential equation [17, 18] with their phase chosen such that the resulting Wannier 
states are real and have the proper parity corresponding to their respective vibrational 
levels [13]. The coefficients a ntq (k) are numerically obtained from algorithms for the 
computation of Mathieu coefficients [19]. Inserting (6) in (3) and taking into account 
the parity of the Wannier states, or equivalently the parity of the band n, one eventually 
arrives at the following expression for the Franck-Condon factors 



27T 

(k + —q)(Ax + r 



a n, q (k)a n , )q (k), (7) 



'(Ax) = 2 J2 5> 

fceBZ q&L 

where we have defined J-"(x) := cos(x) if n and n' have the same parity, and 
J-'ix) := sin(x) otherwise. Numerical evaluation of (7) is shown in figure (2). 

Considering a single lattice site and assuming the harmonic approximation for the 
potential, the shift operator takes the simple form T& x = exp[r] x (a Ji — a)], where a) (a) 
is the raising (lowering) operator acting on the vibrational states. Here we introduced 
the spatial Lamb-Dicke parameter 

Vx = Ax/(2x ) , (8) 

where xq is equal to the rms width of the motional ground state. When r) x <C 1, 
taking the first order term in rj x of T& x allows for a simple expression of the Franck- 
Condon factors for transitions on the same lattice site (i.e., r=r'=0), /™q°(Ax) ~ 
On,n' + r) x {\fnJ8 n > tn+1 - \A^n',n-l)- 
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Figure 3: State-dependent optical lattices relatively shifted by a distance Ax. The total 
trap depth difference AU tot = U^ ot — U^ ot , and lattice contrast W s for spin state \s) are 
shown. Unlike the spin |f) lattice, the contrast and total depth of the spin \\) lattice 
vary with the shift distance. 

2.2. Experimental setup 

We load Cesium ( 133 Cs) atoms from a magneto optical trap into a ID optical lattice 
formed by two counter-propagating, far-detuned, linearly polarized laser beams. The 
filling factor is at most one atom per lattice site due to light-induced collisions [20]. 
A weak guiding magnetic field of 3 G oriented along the lattice lifts the degeneracy 
between the Zeeman sublevels of the Cesium 6 2 Si/2 ground state such that atoms can 
be initialized by optical pumping beams into the hyperfine state \\) = \F = 4, mp = 4). 
Microwave radiation, at around Uuw = 2n x 9.2 GHz, couples states |^) and \\) = 
\F = 3, mp = 3) with the bare Rabi frequency of Q = 2ir x 60 kHz [21]. The spin state 
of the atom is probed using the so called "push-out" technique [22] which consists of 
counting the fraction of atoms left in \\) after all the atoms in |f) have been removed 
by an intense radiation pulse. 

An angle 9 between the linear polarization vectors of the two beams forming the 
lattice is equivalent in the circular basis to a phase delay of 29 between two collinear 
and independent circularly-polarized standing waves, <j + and a~ , or equivalently to a 
standing wave longitudinal relative shift of 



The polarization angle 9 is controlled by an electro-optical modulator (EOM) and two 
quarter-wave plates in the path of one of the two lattice beams. The two in-phase 
circular components of the beam are mapped by the first Al/4 plate onto orthogonal 
linear polarizations parallel to the EOM axes. The retardation 29 induced by the EOM 
is proportional to the voltage signal applied to it. The last plate then converts the linear 
polarizations back into the circular ones while conserving the delay. 

The trapping potentials resulting from the a + and a~ standing waves for a spin 
state \s) are 



Ax sw (9) = 9d/n. 



(9) 



U s = U t s ot + W s cos 2 [k L (x-x° s )] 



(10) 
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where W s and f/* ot are the lattice contrast, taking positive values, and total trap depth 
for state \s), respectively. Both W s and f/* ot depend on the lattice lasers wavelength 
Al and on the lattice shift Ax or equivalently the polarization angle 9, see figure (3). 
For alkali atoms, one can define the "magic wavelength" as the one where the state |f) 
experiences the a + standing wave only. This occurs at Al = A2 + (Ai — A2)/(2Ai/A2 + l) ~ 
A2 + (Ai — A2V3, where Ai (A2) is the wavelength of the Di (D2) line [23, 24], which is 
Al = 866 nm in our case. At this wavelength, for the spin |^) state, equation (10) reads 

U t = -W t + W t cos 2 (k L x - 6/2), (11) 

while the spin \\) state experiences both a + and a~ standing waves with a relative 
weight of 1/8 and 7/8, respectively. The lattice potential in this case is 

U i = -W t + (l/8)W t cos 2 (k L x -9/2) + (7/8) W t cos 2 (k L x + 6/2) . (12) 

With the notation of equation (10), one finds that = —Ul ot is independent from 
the angle 9, while W i = [cos(fl) 2 + (3/4) 2 sin(fl) 2 ] 1 / 2 W t and f/] ot = -(W t + W±)/2. In 
addition, one obtains the lattice relative shift Ax = (d/n){9 + arctan[3tan(#)/4]}/2. 
Equations (11) and (12) constitute the closest realization of the idealized spin-dependent 
lattice discussed in section 2.1. The small admixture of a a + component in equation (12) 
results in a lattice depth W± that depends on 6, or equivalently on the lattice shift Ax, 
which makes the energy levels e Stn (Ax) depend on the spin state and on the shift Ax, 
see figure (3). The nonlinear position shift of the potential, x°, makes Ax deviate 
from the standing wave relative shift Ax sw in equation (9), and this has to be taken into 
account in the calculation of the Franck-Condon factors [23]. 

The typical total lattice depth used in our experiment is ~ 850 E^ tt 

(corresponding to 80 uK), with E^ tt = h 2 k 2 j /2rriQ S as the lattice recoil, which amounts to 
an oscillation frequency along the lattice axis of u v ib ~ 2tt x 116 kHz. In the transverse 
direction, atoms are confined only by the Gaussian profile of the lattice lasers which 
results in a transverse oscillation frequency of u ra d ~ 2rr x 1kHz. The typical initial 
temperature of the atoms loaded into the lattice is T w 10 uK, which in the harmonic 
approximation amounts to mean vibrational numbers of n v ;b ~ 1.4 and n ra d ~ 280 in 
the axial and transverse directions, respectively. 

2.3. Microwave sideband spectra 

We investigate sideband transitions by recording microwave spectra for different lattice 
shifts. Controlling the relative distance Ax allows us to continuously tune the parameter 
f] x from to about 5. In order to resolve the sidebands we use Gaussian microwave pulses 
with a FWHM of 30 us and a bare Rabi frequency of Qo/2it = 36 kHz, corresponding to 
the 7r-pulse condition for the carrier transition. Figure (4) shows a combined spectrum 
where transitions from n = to levels up to n! = 14 are well resolved [11]. Four 
spectra are recorded for four different lattice shifts. With an unshifted lattice only the 
carrier transition is visible, and it defines the zero of the microwave detuning Smw- 
The remaining three lattice shifts were chosen such that for each shift distance Ax the 
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sideband coupling strength on the same site (i.e., r=r'), Q,™q(Ax), is simultaneously 
close to maximum for a small group of adjacent sideband transitions. The coupling 
strength for sites r ^ r' can be neglected at the given shifts. 

For each shift distance Ax the microwave spectra are fitted using the spectra yielded 
by a numerical calculation of the time evolution based on the following Hamiltonian 

H — H + H M w (13) 

with 

H = ^ ^ (e Stn (Ax) + £,,t frvns) \s, n, r) (s, n,r\ , (14) 
s=t4 n 

h M w= -^oEE / ^( Aa; )( e " iW it^^)a^> , i+ h - c -)' ( 15 ) 

r,r' n,n' 




Lattice shift Ax/d Lattice shift Ax/d 



Figure 4: (a) Microwave spectrum of sideband transitions \f,n = 0) -H- \l, n') for lattice 
shifts Ax ={0nm (•), 43nm (o), lllnm(H), 176nm (□)} corresponding to the parameter 
r] x = {0,1.2,3.1,4.9} defined in section 3.1 (data points from [11]). The microwave 
detuning is given with respect to the carrier transition frequency. Data points are the 
average on about 100 atoms and they are here fitted with a model that takes into 
account broadening mechanisms detailed in the text. The error bars, reported only for 
three representative peaks, are obtained with the 68% Clopper-Pearson interval method 
for binomial statistics. The panels (b) and (c) compare the expected values (dashed 
lines) for the lattice contrast and total trap depth difference AU tot (see figure (3) 
and text in section 2.2) with the values extracted from the fits (1% uncertainty). 
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Radial position p (|jm) 

Microwave detuning 



Figure 5: (a) Inhomogeneous broadening effect due to the transverse motion of the 
atoms in the trap. The overall peak profile (bottom curve) is the convolution of the 
other two profiles. The Fourier-limited FWHM and the thermal broadening are typically 
20kHz and 5kHz (for T 2 d ~ 10 pK and first sideband), respectively, (b) Left panel: 
Franck-Condon factors as a function of the radial distance p (here, 9 = 15°). The gray 
profile shows the 2D radial distribution from equation (16) for the same temperature. 
Right panel: Resulting thermal distribution of Franck-Condon factors. 



where c^hs denotes the hyperfine splitting frequency of the ground state. With this 
notation, the microwave detuning reads Smw = wmw — ^hs- 

Given the deep lattice regime considered here, in the numerical solution of equation 
(13) the maximum number of vibrational levels per site can be restricted with a good 
approximation to n max = 15, before atoms start to behave like free particles tunneling 
between sites or directly coupling to the continuum. In this regime, the coupling strength 
for a sideband transition between two lattice sites separated by a distance x > d are 
two orders of magnitude lower than the typical time scales of our experiment; therefore, 
we limit the site indices to r = r' . 

In the fitting of the sideband spectra, the energy levels e s ^ n and Franck-Condon 
factors J"q° depend on the fitting parameters Ax, t/* ot and W s . In particular, in the 
harmonic approximation the spacing between two adjacent peaks is equal to the trap 
frequency of the U± potential, which therefore determines the lattice contrast W±; the 
absolute offset of each spectrum is mainly determined by the difference of the total trap 
depths, AU tot = £/| ot — t/| ot , expressed in frequency units. Additionally, an average 
over the thermal motion of the atoms in the transverse direction of the one-dimensional 
optical lattice has to be performed. In fact, the lattice parameters [7* ot and W s depend 
on the transverse position of the atom, and to take this dependence into account 
we assume that during the microwave dynamics an atom has a "frozen" transverse 
position p. This assumption is justified by the slow transverse motion of the atoms, 
w m d/27r « 1 kHz, compared to the lowest bare Rabi frequency used for the microwave 
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pulse, VLq/2tt m 14 kHz. The transverse positions of the atoms are then assumed to be 
distributed according to a two-dimensional Boltzmann distribution, shown in figure (5b) 
and given in the harmonic approximation by 



with T 2 d the transverse temperature. The thermal transverse position distribution 
results in an inhomogeneous distribution of microwave sideband resonance frequencies 
and of Franck-Condon factors, shown qualitatively in figure (5). Both distributions 
are used to weight the calculated spectra, with T 2 t> as an additional fitting parameter. 
The figure shows that thermal broadening effect becomes larger for higher sidebands, 
exhibiting a more pronounced asymmetric peak shapes. This behavior has a clear 
explanation: in the harmonic approximation, for instance, one expects the thermal 
broadening to increase linearly with band index n while the Fourier-limited FWHM 
remains constant. 

The best-fit results for W± and AU tot are shown in figures (4b) and (4c). This 
method allows us to spectroscopically determine the parameters of the spin-dependent 
potentials seen by the atoms with a relative uncertainty of about 1%. The small 
deviations from the expected values (dashed curves in the figure) can be attributed in 
part to measurement uncertainty and to polarization imperfections in the standing wave 
beams, resulting in slightly distorted potentials. For instance, polarization distortion 
can be responsible for the non monotonic behavior of the data points in figure (4c). 
From the fit, we obtain a temperature of T 2 d = (2.7 ± 0.5) uK. Without axial ground 
state cooling, we measure a three-dimensional temperature of 10 uK by means of the 
adiabatic lowering technique [25]. This discrepancy requires further investigations. 

3. Microwave sideband cooling 

The general principle of resolved sideband cooling, depicted in figure (6), relies on 
the repetition of cooling cycles where each cycle starts by a sideband transition 
|f, n) — > \l,n — 1) removing a vibrational energy quantum hw v ib- The cycle is then 
closed by an optical repumping process with a transition to an optically excited state 
|e) followed by a spontaneous decay to the initial spin state. Because of the optical 
repumping, the motional energy of the atom in each cycle increases on average, which 
corresponds to heating. Therefore, in order to achieve cooling the overall energy gained 
by an atom after one cycle must be negative. In general, heating is caused by the 
momentum recoil from the optical repumping photons, i.e. recoil heating. In the 
microwave-based scheme however, shown in figure (6b), an additional source of heating, 
called hereafter "projection heating," is present. It is due to the difference between the 
trapping potentials of the internal states, in this case it is a spatial shift. This difference 
makes the projection of a vibrational state \\,, n) on an arbitrary state |f, m) appreciable 
in contrast to the case of identical potentials where transitions beyond m = n, n ± 1 are 




with 




(16) 
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Figure 6: (a) Raman sideband cooling scheme: a two-photon Raman transition between 
two identical trapping potentials reduces the vibrational state, |t, n) — » \l, n — 1). The 
wavefunction is shifted in momentum space by hAk. (b) Microwave sideband cooling 
scheme: a microwave transition between two shifted trapping potentials reduces the 
vibrational state. Note that we use here a blue sideband transition to reduce the 
vibrational state, instead of the typical usage of a red sideband transition [7]. 



negligible. 

3.1. Raman vs. microwave sideband cooling 

In the standard Raman-based sideband cooling schemes [6, 26] the sideband is induced 
by a two-photon transition where the coupling is given by the matrix element 

fin~ = (l,n- l \T Ak \ t, n) x fi*™ , (17) 
where T^fe = exp(ixAk) is the momentum shift operator and Ak ~ 2k opt is the 
wavevector difference between the two optical photons for counterpropagating beams. 
From here on, it is understood that all transitions occur on the same site, r = r'. In 
the microwave-based scheme we neglect the microwave photon recoil, and the sideband 
coupling corresponding to a lattice shift Ax between nearest neighboring sites is then 
given by 

^n— l,n 

(i,n-l\T Ax \ t,n) xO . (18) 

Using the harmonic approximation, the Raman and microwave sideband couplings can 
be expanded to the first order in the parameters rjk = frk opt /(2p ) and rj x = Ax/(2xq), 
as shown in table (1), where po = ^mc s hjj v ^/2 and xq = \Jh/ (2mc s w v ib) are the 
momentum and spatial rms width of the ground-state wavefunction, respectively [27]. 
From table (1) we can note a clear duality between momentum and spatial shifts in the 
two sideband cooling methods. The duality is better emphasized by using the general 
complex Lamb-Dicke parameter 

r] = r] k + ir] x = hk opt /(2p ) +iAx/(2x ) = k opt x + ip Ax/H, (19) 
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Table 1: Raman vs. microwave sideband cooling in the harmonic approximation. The 
sideband couplings are first-order expansions in rjk and r\ x in the Lamb-Dicke regime 
defined by \r)\ = |% + if] x \ ^ 1> under the harmonic approximation. 





Raman 


Microwave 


Sideband coupling strength fi n _i jn /fio 


i2i] k y/n 


-VxVn 


Recoil heating per cycle 


2hu vih rf k 


2huj vih ril 


Projection heating per cycle 






Overall heating per cycle 


2kj vih T]\ 


^w vib [ril + 2ri 2 k ) 



which accounts for both degrees of freedom via the momentum and spatial Lamb-Dicke 
parameters, rjk and rj x , respectively. This generalized approach was introduced first in 
ion systems to describe microwave-induced sidebands in the presence of spin-dependent 
forces [14]. 

In the Raman-based cooling schemes with identical trapping potentials, the spatial 
Lamb-Dicke parameter ri x vanishes and the heating comes from the recoil of the optical 
repumping photons, as depicted in the figure (6a). In the microwave-based scheme 
however, the generalized Lamb-Dicke parameter is complex and the heating is caused 
by a combination of recoil and projection heating. 

The energy gained by an atom from recoil heating after one cycle results from two 
recoils, one from absorption and one from spontaneous emission, and is therefore given 
by 



AE r . 



2E 



R, 



(20) 



where = h 2 k^ pt /2mc s is the optical photon recoil energy [28]. This quantity does 
not depend on the details of the potentials but only on the atom's properties, and it 
expresses the overall three-dimensional recoil heating. 

In the shifted potentials shown in figure (6b), in addition to the recoil heating, 
the atom's motional energy increases on average by the projection heating energy 
AE pTO y This is due to the non- vanishing projection of the atom's initial vibrational state 
\-l,n) onto the vibrational basis |1\ m ) of the final spin state in the optical repumping 
process. In the harmonic approximation, with H ext = /za; vi t,(n + l/2), and after adiabatic 
elimination of the excited state |e), the projection heating contribution for a relative 
shift Ax can be derived as 

2 



AE- 



proj 



m 



T 



fi^vib S ^{ m - n) 

m 

^2 ( m | [HT Ax - f Ax H 

n 

n I (H ext (Ax) - H, 



Ax 



11 



11 



n 



T 



Ax 



in 



ext ) 



11 



(21) 
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where we have introduced the shifted Hamiltonian 

v 2 1 

H ext (Ax) = r Ax H e ^T Ax = ^— + -m Cs u4 b (z + Ax) 2 . (22) 

zmcs £ 

The result of equation (21) applies in general for any potential profile, and in the 
harmonic approximation it results in a quantity which is independent of n, 

A ^proj = l,mcsU 2 ih Ax 2 , (23) 

which is nothing but the potential energy difference as expected from the semi-classical 
picture in figure (1). 

Using the same method, one can generally show that in the microwave sideband 
cooling scheme the total average heating energy gained by an atom in one cooling cycle 
is the sum of the recoil and projection contributions. The total energy balance per cycle 
then becomes 

AE tot = AE pm j + AE iec - hw vih = huj vih (r] 2 x + 2r] 2 k -l). (24) 

Similarly to the usual definition of the Lamb-Dicke regime [29], the condition for cooling 
AE tot < defines a generalized Lamb-Dicke regime as the range where \rj\ < 1. 

3.2. Quantitative model based on master equation 

The general theory of sideband cooling is very well known and has been extensively 
studied in the literature [5, 27, 30, 31]. Here, we discuss a quantitative model based on 
the Lindblad master equation formalism. To provide a concrete example, we apply the 
model to the level scheme of our specific system, though the model can be adapted to 
other similar systems. 

In the cooling cycle depicted in figure (7), microwave radiation resonant with the 
first blue sideband transfers atoms from states |t, n) to states \\., n — 1). Concurrent 
with the microwave, a <r + -polarized repumper laser beam couples state \\) to state 
|6 2 -P3/2, F' — 4) = |e), from where the atoms close the cooling cycle by spontaneously 
decaying back to state Due to the appreciable probability of atoms decaying from 
state |e) to state \F = 4, mp — 3) = \a), a second equally polarized pumping laser 
couples the two states and brings the atoms which have decayed to state \a) back into 
the cooling cycle. In each cycle, an atom loses energy on average until it reaches the 
"dark state" |j, n — 0) where it is no longer affected by the microwave or the repumping 
lasers. Nevertheless, a small probability remains that the dark state is depopulated 
due to photon scattering from the lattice lasers or an off-resonant microwave carrier 
transition. 

To describe the cooling dynamics, we reduce the problem at hand to an effective 
model with three spin states with the set of motional states associated with each one of 
them. The considered Hilbert space is then the one spanned by the states \s,n), with n 
being the vibrational level and \s) being one of the three internal states |t) ?|4) or l a )- 
The optically excited state |e) is adiabatically eliminated due to its very short lifetime, 
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t = 30 ns, compared to the motional time scale. We use the Lindblad master equation 
formalism to write the time evolution of the effective model's density matrix as [31] 

^= -^[H' + Huw,p]+C[ P ] : (25) 
where H' Q is the extension of the Hamiltonian from equation (14) to the states \a, n) 

^o= Yl Y £ sA Ax ) \s,n,r) (s,n } r\ , (26) 

s={t,-U} n 

and C is the Lindblad superoperator with the projectors 

L<£' = \s,n,r)(s'yy\, (27) 

and the effective decay rates jg'^'f' for the transitions \s',n',r') — > \s,n,r) which are 
given by 

itnf = * s R s ,(\M s : ; :'f\ 2 ) , (28) 

with 

M^'f = (n, r, s\ f Aks s , f Ax \s', n', r') . (29) 

Here, a s is the branching ratio for the spontaneous emission from state |e) to state \s), 
and R a , Ri are the pumping and repumping rate, respectively, as shown in figure (7). In 
addition, we account for the possibility that an atom in |^) state scatters a photon from 
the lattice with the rate R^. The matrix element Mf'^'f' accounts for the relative spatial 
shift between the two involved vibrational states and for the transferred momentum of 
both optical photons A.k StS > = k opt + k sp ■ e x in the optical repumping process, with 




Figure 7: Microwave sideband cooling scheme in a realistic physical system using 
133 Cs atoms, (i) Microwave radiation tuned to the first blue sideband induces a 
|t, n) — > \l,n — 1) transition decreasing the motional quantum number by one. (ii) 
The cooling cycle is closed by an optical repumping transition \\) — > \F' — 4), with 
rate R^, and (iii) a spontaneous decay back to state In (iv) an additional pumping 
laser brings the atoms which have decayed to state \a) back into the cooling cycle, with 
rate R a . Atoms reaching the dark state |f, n = 0) are out of the cooling cycle unless 
off-resonantly excited or heated externally. 
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k sp being the wavevector of the spontaneously emitted photon and e x being the unit 
vector along the lattice direction. Additionally, one has to perform an average over k sp , 
indicated by the angle brackets in equation (28) . 

Given our experimental parameters, we compute the steady-state solution to 
equation (25) numerically, using the same approximations as in section 2.3. In the 
computation, the microwave is resonant with the |^, 1) — > \\,, 0) transition. The rates 
R a and R± are set by the experimental values, which are chosen comparable to fi ,i and 
smaller than the vibrational level separations to avoid off-resonant transitions by power 
broadening of the vibrational levels of the F = 3 ground state. Figure (8) shows a 
contour plot of the ground state population P\ n =o) = r P\s,o,r) i n the steady state as a 
function of the bare microwave Rabi frequency and the relative shift distance expressed 
in terms of rj x . When projection heating dominates, rj x > rj^, the energy balance in 
equation (24) just requires r] x < 1 for cooling; for instance, figure (8) shows that a 
ground state population P| ) > 80% can be reached with r\ x < 0.8. For very small lattice 
shifts however, with r) x <C 1, the microwave coupling for the blue sideband transition 
becomes small compared to that of the carrier transition, which renders the microwave 
action of removing an energy hu) v fo per cycle inefficient compared to the recoil heating, 
which is the dominant heating source for such small shifts. Weak microwave sideband 
coupling and hence inefficient microwave cooling will also be present at very low bare 
Rabi frequency, namely at the Rabi frequencies where the sideband coupling becomes 
lower than the rate of depopulation of the dark state. For high Rabi frequencies of the 
same order of magnitude as the vibrational level spacing, the carrier coupling becomes 
comparable to the blue sideband coupling, and the microwave cooling action is again 
reduced. 




Vx 

Figure 8: Steady state population in the motional ground state P\ n= o) as a function of 
r] x and the bare microwave Rabi frequency Qq. 
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3.3. Experimental results 

Microwave cooling is obtained by applying microwave radiation on resonance with the 
first blue sideband, |f, 1) — > \\,, 0), for a certain duration r coo ii ng , at a certain lattice shift 
Ax, concurrently with the two optical pumping lasers as shown in figure (7). In order 
to probe the final vibrational state distribution, we record a spectrum of the first order 
sideband transitions using a Gaussian microwave pulse satisfying the 7r-pulse condition 
for the first red sideband, figure (9a). In the low temperature limit, the height of the first 
blue sideband peak provides a good measure of the motional ground state population, 
i"W o) ) an d thus of the cooling efficiency For instance, for atoms in the ground state one 
expects to detect no blue sideband. Figures (9a) and (9b) show two microwave spectra 
recorded before and after cooling, clearly indicating a reduction of temperature by the 
cooling process. 

In order to determine the optimum cooling parameters, the blue sideband height is 
remeasured while scanning different variables, namely the optical pumping intensities, 
the cooling microwave power and frequency, the lattice shift distance and the duration 
of the cooling pulse. Figure (9c) shows a scan of the cooling microwave frequency. As 
indicated by a nearly zero detected signal from the blue sideband, the optimum frequency 
for cooling lies evidently in the vicinity of the first blue sideband, while a less pronounced 
cooling is also present at the position of the second blue sideband. Furthermore, the 
measurement reveals the absence of the blue sideband signal in a broad range extending 
to negative detunings in addition to a weak dip at the position of the carrier. These two 
observations are correlated with a decrease in the atom survival given in the same figure. 
This shows that, instead of being due to cooling, the absence of the blue sideband here 
is due to increased atom losses. In fact, for zero and negative microwave detunings, 
that is if the microwave is resonant with the carrier \t,n) — > \\., n) or red sideband 
|t, n) — > |4, n + m) transitions respectively, the energy of the atom increases on average 
in each cooling cycle. In the case of zero detuning the increase is due mainly to recoil 
and projection heating in the absence of microwave cooling, while for negative detunings 
microwave sideband heating occurs in addition to the recoil and projection heating. 

Once the optimum cooling parameters have been determined, we extract the 
achieved steady state temperature assuming a thermal Boltzmann distribution and 
neglecting the anharmonic spacing of the vibrational levels. The ratio between the 
red and blue sideband heights is proportional to the Boltzmann factor which is related 
to the average motional quantum number (n) by 



Using the fitted sideband heights from figure (9b), we calculate (n) = 0.03 ± 0.01, and 
a ground state population of P^ — 97%, corresponding to a temperature T ~ 1.6 uK. 
Table (2) summarizes the optimum cooling parameters for our setup. 




(30) 
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Table 2: Optimal microwave cooling parameters for w v ib = 2ir x 116 kHz. 
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4. Motional state control 

^.1. Motional state detection 

We have developed a vibrational state detection scheme which allows us to determine 
the vibrational state distribution of any given motional state. It relies on removing all 
atoms above a selected vibrational state n from the trap and counting the remaining 
atoms, as illustrated in figure (10a). The distribution is then reconstructed from the 
differences of subsequent measurements. 

Atoms are first transferred to state |^) by means of an adiabatic passage microwave 
pulse that is resonant with the carrier transition in unshifted lattices, which preserves 
vibrational states' populations. A microwave pulse resonant with the red sideband 
\i,n) — > |1\0) transfers atoms from states m) with m ^ n to states | 7,771 — n). 
The transferred atoms are eventually pushed out of the trap (see section 2.2). However, 
since the sideband transition rates depend on the initial vibrational state \l,n) (due 
to, e.g., trap anharmonicity and Franck-Condon factor differences) the microwave pulse 
does not achieve full transfer efficiency for all transitions. To overcome this limitation, 
the procedure of microwave pulse plus push-out is repeated several times to deplete all 
vibrational states \l, m) with m ^ n. If / is the population transfer efficiency for a 
given n, then after iV repetitions the effective population transfer efficiency becomes 
/' = 1 — (1 — f) N . For instance, an initial efficiency of / = 70% is thus increased to 
/' ~ 97% with N = 3 repetitions. Measuring the fraction of remaining atoms as a 
function of the microwave frequency, we obtain a sequence of plateaus at the successive 
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Figure 9: Microwave spectroscopy performed (a) before cooling and (b) after 20 ms 
of microwave cooling, with optimal experimental parameters (see Table 2). (c) Blue 
sideband height vs. detuning of the microwave cooling frequency (o) and atom survival 
probability measured after the cooling (□), (data points in (a) and (b) are from [11]). 
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sideband frequencies, as shown in figure (10b). The plateau corresponding to the n 
sideband indicates the integrated population of states m < n, that is, the cumulative 
distribution function F n = Y^=aPmi from which the individual populations of the 
vibrational states are then derived. 

Motional state n 




Figure 10: (a) Method for measuring vibrational state population distributions: (i) an 
initial microwave pulse resonant with the n th red sideband transfers all atoms in states 
IJ,, m), with m ^ n, to state (ii) the transferred atoms are pushed out of the lattice; 
(i) and (ii) are repeated N times to overcome low pulse efficiency, (b) Surviving fraction 
of atoms for a thermal state, with the dotted lines indicating a thermal distribution of 
T w 11.6 uK; this temperature is compatible with the independently measured one of 
about 10 uK. For each sideband n, after N = 3 repetitions of the microwave pulse plus 
push-out, only the atoms in states m < n survive. The horizontal dashed line indicates 
the maximum survival probability, limited by the off-resonant transitions during the 
repeated pulses. For the sake of clarity, error bars have been displayed for the carrier 
transition only. 



4-2. Motional state engineering 

With 97% of the atoms cooled to state |t, n = 0) (see section 3.3), controlled preparation 
of different motional states is possible using a combination of microwave pulses and 
selected lattice shifts. 

The simplest state that can be prepared is the Fock state |4,, m). It requires 
addressing the m-th red sideband transition at the lattice shift Ax chosen to maximize 
the coupling |t, 0) -B- |4,, m). The fidelity for preparing this state is limited by the cooling 
efficiency, the population transfer efficiency and the selectivity of the microwave pulse. 
Using an adiabatic passage pulse [32], a state preparation fidelity close to 98% has been 
achieved for states up to m = 6. 

A superposition of two Fock states is created by a two-pulse sequence as shown 
in the inset of figure (lib). An initial microwave pulse resonant with the transition 
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(a) Superposition state (b) Coherent state | a) 




Vibrational state n Vibrational state n 

Figure 11: Motional state preparation and analysis. Shown are the populations of the 
vibrational states n — 0, .., 3 after (a) creating superposition states of |n = 0) and \n = 2) 
with different weights (from top to bottom, area of the first MW pulse 0.30, 0.40, 0.55, 
0.70 in units of n) and (b) coherent vibrational states for different amplitudes a, where 
the left bars (brighter red) indicate the theoretically expected populations. The analysis 
technique used here, see figure (10a), can only measure vibrational states' populations 
but not coherences. 



|t, 0} — > \-l,2), performed at the lattice shift which maximizes the coupling for the 
transition, generates the state 

|V) = ct,o|t,0)+c w ||,2) (31) 

with variable coefficients c^o and determined by the pulse duration. The lattice 
shift Ax is then changed to the distance at which the Franck-Condon factor for the 
transition |t, 2) -H- \\,, 2) is zero. The shifting is precisely timed so that the probability 
of changing the vibrational state by the acceleration of the lattices is zero [33]. At the 
new lattice shift, a second microwave pulse resonant with the carrier transition maps 
the population \c^ } o\ 2 onto \i,0). One would expect as a result a coherent superposition 
between \\., 0) and \\., 2). However, because of the appreciable duration of the sequence of 
320 us (two sideband-resolved pulses plus lattice shift operation) compared to the total 
spin coherence time of ~ 250 us in our setup, the coherence between the two vibrational 
states is partially lost during the preparation procedure. This is a technical limitation 
which can be overcome by improving the coherence time, for instance in our setup, by 
cooling the transverse motion of the atoms to the three-dimensional ground state [34]. 
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This scheme represents a relevant step towards the use of the vibrational state as the 
physical carrier for a qubit and/or the preparation of arbitrary motional superposition 
states when working with neutral atoms. 

In the same vein of engineering motional states, we project the state |t, n = 0) onto 
a shifted potential to create the state 

|«) = f Ax |t, n = 0} = e aat ~ a * a |t, n = 0) (32) 

with a = r] x . We realize this by applying an optical repumping pulse while the lattice 
is displaced by Ax. This corresponds to exciting the transition \\) — > |e) followed 
by a spontaneous decay to state |t), which occur on a time scale much shorter than 
the oscillation period of the atom in the trap. We also neglect the recoil transferred 
by the optical repumping photons, which is equivalent to assuming r)k = 0. Because 
the decay process additionally involves transitions to states \a) and |t), the resulting 
state is a statistical mixture of the three internal states; our analysis scheme however 
measures the vibrational population of the state projected on |t) state shown in equation 
(32). The statistical mixture can be avoided by replacing the optical repumping 
pulse and spontaneous decay by a fast two-photon Raman transition. Measuring the 
population distribution of the created state reveals a clear agreement with the theoretical 
expectation, as shown in figure (lib). With the state detection scheme presented in 
section 4.1, so far we can only measure populations, while coherences could be accessed 
in the future through interferometric schemes. 

5. Conclusions and outlook 

We have shown that microwave sideband transitions in spin-dependent optical lattices 
are a favorable alternative to Raman transitions for sideband cooling and motional 
state engineering. The effective Lamb-Dicke parameter can be continuously adjusted 
from zero to above one, giving the possibility to address directly higher-order sideband 
with coupling frequencies comparable to the bare Rabi frequency. We investigated the 
performance of microwave sideband cooling in the generalized Lamb-Dicke regime, and 
we compared it to the Raman sideband cooling; our analysis can be easily extended to 
the three-dimensional case [15]. 

Quantum engineering of motional states represents one of the most attractive uses of 
microwave-induced sidebands. We demonstrated here a first step towards the creation 
of superposition between Fock states, and the preparation of coherent states. In the 
future, the interest resides in proving the coherence properties of these states through 
interferometric schemes, for instance, by measuring the accumulated phase between two 
distinct Fock states, or through quantum beat experiments [35]. 

Along the same line, spin-dependent shift operations can be employed to transfer a 
state-dependent momentum kick, allowing the realization of a superposition of opposite 
coherent states, producing Schrodinger-cat-like states as has been realized with ions [36] . 

Microwave control of atomic motion in a spin- dependent optical lattice can be of 
interest for storing and processing quantum information via the motional states [37]. For 
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instance, the strength of coherent collisions for atoms close to the motional ground state 
exhibits a marked dependence on the relative motional state, which can be exploited, in 
analogy to [24] , to realize maximally entangled states in the motional degree of freedom. 

In addition, microwave sideband transitions open the way for quantum transport 
experiments, where continuous tunneling between adjacent lattice sites occurs when Ax 
is close to d/2, i.e., close to half the lattice spacing [38]. 

Finally, it is worth noting that the microwave cooling technique studied here does 
not strictly require the use of the "magic" wavelength for the lattice potential, but can 
still be operated with the same efficiency at other wavelengths, e.g., at Al = 1064 nm 
as we have tested. In fact, the optimal cooling efficiency occurring at around t] x ~ 0.3, 
see figure (8), can be reached by adjusting the polarization angle 9. 
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